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1 Introduction 



Although the heavy fermion (HF) behavior of the rare earth- and actinide-based intermetanides 
is related to the localized electrons in the 4f- and 5f-shells of Ce and U, the influence of the 
crystal field (CF) splitting of the f-levels is inessential, as a rule, for the low-energy excitations 
in these systems, since the inequality A^^ T* is valid for most of these compounds [|^. 
Here T* ~ Tk is the temperature which characterizes the energy scale of the HF spectrum. 
The only role of this interaction is to form the ground state Kramers doublet which is the 
eventual source of the low-energy excitations with extremely high density of states. However, 
the picture changes radically when this inequality is violated. One can expect that in case 
of IS.CF < T* the inevitable interplay between the local CF excitations and itinerant HFs 
should result in appearance of the fine structure in the low-energy spectrum which should be 
characterized by at least two additional energy parameters, i.e., the position of the CF level 
and the magnitude of intermixing. This interplay should result in remarkable changes in the 
spectral and thermodynamic properties of the system in the temperature interval < T < Aq f , 
although it is clear that the characteristic HF behavior of the system should persist at T — > 
0: the Fermi liquid-like temperature dependences for the specific heat, C(T) ~ T, electrical 
resistivity, p(T) — p(0) ^ T^, NMR relaxation rate, TiT= const, etc., should be observed, 
although the coefficients in these laws should be sensitive not only to T* but also to these new 
parameters. On the other hand, the hybridization of the HF and CF excitation can result, 
as it frequently occurs, in appearance of the psedogap in the energy spectrum, although this 
pseudogap should be characterized mainly by the crystal field and mixing parameters. 

The idea of applying these considerations to description of unusual properties of CeNiSn 
was proposed in Q . It was shown that one of the necessary preconditions for realization of this 
mechanism is the specific picture of CF level splitting: it was supposed that two lowest Kramers 
doublet of Ce(f^) configuration are separated by a small energy gap 

AcF < T* (1) 

while the distance between the ground state and the second excited level is much larger, ^ 
AcF- The characteristic HF energy T* is estimated as 40 to 80 K from various experimental 
data for CeNiSn Later on, the indirect measurements by Alekseev et al 0] confirmed the 
assumption (Q). In that paper the crystal field was measured on the Nd ion in Ndo.3Lao.7Ni 
crystal which is equivalent to CeNiSn from the point of view of RE crystalline environment 
and lattice spacings. The crystal field acting at the RE site was restored from measured CF 
splitting of Nd f-level and then recalculated for the case of Ce ion with the same parameters 
of equivalent CF Hamiltonian in approximation of purely electrostatic field. The CF splitting 
energies were estimated as Aqf ~ 4.4 meV and A^^ ~ 14.0 meV which is consistent with the 
assumption of ref . Q and the inequality ([l|) . 

Another experimental result crucial for the above general consideration was the observation 
of Fermi-liquid behavior at low enough temperature in good quality samples of CeNiSn and 
CeRhSb. It was found that these samples demonstrate the metallic T^-dependence of electrical 
resistivity at T < 7K |5|, the NMR relaxation rate returns to Korringa law at T < IK and the 
limiting value of the Sommerfeld coefficient 7 in the linear-T term of C (T) is 7 = 40 mJ/K^mol 
at 0.03 K m. This value is noticeably lower than that predicted by the high-temperature 
estimations of T* . These results seemed to be striking, since the semiconductor-like rise of the 
resistivity and dependence of 1/Ti was observed earlier on the less perfect samples [Q]. 
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The most challenging experimental observations on CeNiSn are the inelastic magnetic neu- 
tron scattering spectra which demonstrate extremely complicated (Q, cj) -dependent structures 
with quasi gap behavior for some directions of Q ^, |l^. All these features disappear 

completely at T=25 K. This picture is inexplicable in terms of standard notions of intersite spin 
fluctuations. At the same time, these data are highly informative from the point of view of the 
structure of low energy spin excitations, so the adequate interpretation of the neutron scattering 
data is the most effective criterion for selecting out of various theoretical interpretations. 

In the present paper the spectral properties of CeNiSn are analyzed on the base of the ideas 
formulated in Q and applied for the quantitative explanation of low temperature thermody- 
namics in ||2|, |lT[. We apply these ideas to the spin liquid which is formed in Kondo lattices, 
provided the magnetic order does not appear at low temperature. We describe this spin liquid in 
terms of resonating valence bond (RVB) excitations similar to those introduced by P. Anderson 
and used later in description of nearly magnetic state in 2D Cu-0 planes of high-Tc materials 
|12|. In this case the characteristic properties of HF liquid are determined by the subsystem 
of spin excitations (spinous) which obey Fermi statistics at low temperature. The stabilization 
mechanism for the RVB state in 3D Kondo lattices was offered in |13|. We modify the theory 
of spin liquid for the case of IS.cf < Tk, calculate the spectrum of RVB excitations for a spe- 
cial case of orthorhombic CeNiSn lattice and give a detailed comparison of calculated neutron 
scattering cross section with the experimental results of P, 10|. The quasi 2D character of spin 



excitations discovered in these experiments is essentially used in our calculations. Quite rea- 
sonable agreement between the theoretical curves and experimental constant-Q and constant-ii^ 
scans seems encouraging. To make an additional test of adequacy of the theory, we calculate 
also the low-temperature specific heat of CeNiSn within the same set of model parameters. 



2 Degenerate spin liquid in a crystal field 

The theory of HF systems starts with the Anderson lattice Hamiltonian. In a special case of 
Ce(^+)(/^) ions in a crystal field this Hamiltonian is usually written as 

H = Hf + H, + H,f + H' (2) 

Here Hf describes the Ce ions in a crystal field, 

^/ = E ^r//A/iA + y E Zilv/iA/ilv'/iA' , (3) 
i,A i.A^A' 

/j'^ is the creation operator of the f-electron in a state A = Tv in a site i, where v is the row 
of the irreducible representation F of the crystal point group, Ey denotes the energy of the 
f-level in the electrostatic crystal field, and U stands for the Coulomb repulsion between the 
f-electrons in the same site. The f-ions are immersed in a Fermi-sea of conduction electrons 
which are described by the Hamiltonian 

Hc = ^ekvclf^CkA. (4) 
fc,A 

The third term in the Hamiltonian (^) describes the hybridization between f- and conduction 
electrons, 

= E E E (^/iACfcA + H.c) (5) 

kk' i A 
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The so called Coqblin-Cornut (CC) approximation ^] which represents the Bloch functions 
by their partial waves c^^, and takes into account only the diagonal in A hybridization matrix 
elements V^^ = (kA|l/|iA) is used. All the one-electron processes beyond this approximation 
are collected in the last term H' in the Hamiltonian (|2|). 

It is widely believed that the s=l/2 approximation for the total moment j of Ce(f^) ion works 
well in the Anderson lattice Hamiltonian because the CF splitting of the sextet j=5/2 usually 
results in a simple Kramers doublet ground state, and the CF splitting exceeds essentially the 
characteristic interaction energy Tk- This procedure seems to be reasonable when one deals 
with the ground state and the low-energy excitations with characteristic frequencies Tiui < Acf- 
Since our task is inclusion of CF excitations in the general picture of spin liquid behavior, we 
should begin with generalization of s=l/2 theory for the case of realistic crystal field scheme. 

The Kondo lattice deals with the well localized f-electrons for which the inequality VkA ^ 
ep — Ef, where Ej is the energy level of free Ce(f^) ion, is assumed to be valid, and the 
Hubbard parameter U is taken large enough to suppress the doubly occupied f-states of this 
ion. In this case the canonical transformation eliminating the hybridization interaction from the 
Hamiltonian (|2|) can be done (see fl^), and the effective Hamiltonian projected to the subspace 
of homopolar states with nj = 1 can be written as 

H^cc) ^Hf + H, + He, + Hh (6) 

where 

^e. = E E E jt'^'iK k')flU,4^,CuK (7) 
kk' i A,A' 

describes the effective sf-exchange interaction, and 

Hh = -Y.Y.Y.Jt''{k,k)flU. (8) 

k i A 

corresponds to effective covalent contribution to the one-site CF splitting due to virtual sf- 
transitions. Here 

jAA'(fc,fcO= ■ (9) 

e/c - Ef 

We assume for the sake of simplicity that the conduction electrons are degenerate in A near the 
electron Fermi surface and neglect the CF splitting of the f-level in the denominator of effective 
sf-exchange integral. 

We are interested in the case of n/ = 1 when the charge fluctuations in the f-channel are 



completely suppressed. Then, according to the scenario of ||13|, |16|] the neutral spin liquid should 
arise at some temperature T* ~ Tk instead of antiferromagnetic or Kondo-singlet state offered 
by the standard Doniach's dichotomy. And, as a result, the spin liquid excitations together with 
the low-energy electrons turn out to be responsible for the low-temperature thermodynamics of 
the Kondo lattice. Now we incorporate the CF excitations in this picture. The "fast" electrons 
with characteristic energies hco > Tk are responsible for formation of the spin-spin correlation 
at T > Tk (the Kondo temperature is defined as Tk ~ e/ exp(— l/2a), where a = JgfMi^ep), 
Jgf = \Vkp\'^/{eF — Ef), and M{eF) is the electron density of states at Fermi level). These fast 
electrons can be integrated out, and the effective Hamiltonian for the spin degrees of freedom 
can be written in the following form. 

= Hf + H'f^ + H^RKj^Y + Hrkky ' (10) 
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i^i' 

HrkKY = H X] ^S'^fiGufiGu'flcu'fi'Gu , (H) 
ii' uu' 

i^i' 

HrKKY = H H I?/ f-iGuf-^Ev' fl Ev' h'Gv (12) 
ii' Ei/!/' 

Here the notations A = Gv and A = Ev are used for the ground state and excited states, 
respectively, H'^ includes the scattering correction to the CF level renormalization due to the 
interaction H^.^, 

4^ ~ a2^G(kF,Ri - ni')BG{9,4>)K{T) (13) 

is the indirect exchange interaction which contains oscillating RKKY function ^gO^f, R-i — R-i')' 
anisotropy factor Bg{9,(I)) iO is the angle between the z-axis which is assumed to be the axis 
of magnetic quantization, and the line connecting the ions and Rj' and (j) is the angle of 
rotation about z-axis) (see Appendix), and additional enhancement factor K{T) due to the 
high-temperature one-site Kondo scattering |13]. I^,^ is the integral of the same type as that 



in eq. (|T^), with its own factors Bqe and ^ge, and Kge{T) = 1. 

The uniform spin liquid state of RVB type in a standard s=l/2 Hesenberg model is described 
in terms of correlators [^], 

A = E(/il/iv), ^/i'- (14) 

It was shown in |]l^] that the weak-coupling Kondo interaction at T > Tk favors the stabi- 
lization of RVB state against antiferromagnetic ordering because the Kondo scattering screens 



the local moment leaving intact the singlet RVB correlators (14). This stabilization "quenches" 
the Kondo scattering at temperatures T ~ T* > Tk and allows the description of spin liquid 
in terms of RVB correlation functions at T < T^-, since Tk is no more a singular point of 
perturbation theory. Having in mind this stabilization mechanism we use the mean-field ap- 
proximation for the uniform RVB state even at T < Tk- Although the interaction with "slow" 
conduction electrons and the influence of spin and gauge fluctuations can modify the properties 
of spin-fermions, we use this mean-field description as the basic approximation for studying the 
system with the CF excitations involved. 

Thus, the correlator similar to (|l^) can be introduced for the ground state doublet Gv: 

A'' = Y.{flaJi'Gu),iy^i'- (15) 
The spin-fermion excitations which are described by the Hamiltonian (10) are constrained by 



the condition 

E/iA/iA = l- (16) 

A 

Important contributions to the interplay between the low-energy spinon excitations and the 
one-site CF excitations are connected with the interactions B' which appear in the Hamiltonian 
(III) beyond CC approximation. The origin of these interactions is the hybridization V^'^ = 
(iA|y|A;A') where V is the component of the crystal field which has the symmetry lower than 
that diagonalizing the energy terms Ey- When the integrals V^'^ are included in the canonical 
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transformation, additional effective exchange terms appear: 

t7i-B*T/i T/i* t/iG 

^kG ^k'G f] f ^+ ^. I ^kG^k'E A f ^+ ^ j_ u 
^ JiEunGu'C,^,Q^,CkGu" + ^JiGuhGu'Cf,,j^^„CkGu + ^La 

\n) 

Tliese terms being combined with CC exchange term ^ gives additional contributions to the 
intersite effective exchange, and the first order terms in V have the following form, 



HtKKY = T. E [lv/fL'hE.flGuh'G."+H.c] . (18) 



Here I^,^ is the integral of the same type as that in eq. (|T3|), but the corresponding coupling 
constant a' is proportional to V^V and the factors Bqe and K{T) differ from those in (jl^). 

Then, introducing the correlator A*^ (|l5|), we obtain the mean- field Hamiltonian for the 
spin-fermion spectrum [|: 

i^i' 

Hmf = EgY. fiGuhcu + E E '^i'/L/i'G. + 5Hmf , (19) 

i,v ii' f 

where 

5Hmf = AcF E fLfiEu + E [^S^fLfiEW + H.C.] +Y,Y. [B^^fhuhEW + H.c] (20) 

'\u iuu' ii' i/t'' 

We use the bare ground state CF level as a reference point, Eq = 0, then 

Eg = -N~' jF''iKk) + Y.4''^B^ + Bg (21) 

k,k>kp i' 

includes the covalent and indirect exchange contribution to CF shift determined by the inter- 
actions dH) and ([lT|) , respectively, and 

AcF = AcF-N-' ^ Jf^(k,k)^i?f (22) 

also includes the one-site covalent corrections. The coupling constant in the MF Hamiltonian 
is given by the following equation, 

7,v = T.lS'''{fLhc.)=4^A''. (23) 

Only the sites connected by antiferromagnetic RKKY coupling contribute to the RVB correla- 
tion function. 

The one-site and intersite mixing constants are determined as 

fig^ = -V (24) 



^The non-diagonal term ( p^ gives no contribution to the uniform mean-field RVB pairing and describes only 
the fluctuation corrections to the mean-fleld solutions. These corrections are beyond the framework of the present 
paper 



6 



and 

= It/Af" (25) 

respectively. 

Then, one should diagonalize the quadratic form ([T9|), ( pO[ ) and obtain the spinon dispersion 
law 6/3 (k) {(3 are the indices of spinon band) under the global constraint condition 

(2iV)-ig(l-tanh^^^^)=l (26) 

which is used in the mean- field approximation instead of exact local constraint ([T^). This 
equation determines the spinon chemical potential /x. 

Thus, to find the spin-fermion spectrum in the mean-field approximation we need only to 
diagonalize the Hamiltonian (p!9|), (pO|) under the constraint (|26 



3 Spinon spectrum of CeNiSn 

CeNiSn crystallizes in the orthorhombic lattice which belongs to the noncentrosymmetric space 



group Pn2ia |17]. This structure can be described as zigzag chains of Ce atoms directed along 
a axis (easy magnetization axis) and surrounded by slightly distorted trigonal "drums" formed 
by Ni and Sn ions. Thus, the point symmetry of the crystal field on Ce ions can be treated 
as nearly trigonal (-Dad) with rotation axis parallel to a-axis of the crystal, and the monoclinic 
distortion {Cg) can be considered as a small correction to the trigonal crystal field Q. Each 
elementary cell contains four Ce ions. 

To apply the theory of spin liquid state explicated in Section 2 to the case of CeNiSn one may 
use the irreducible representation of the trigonal point group D^d as a basis for Cornut-Coqblin 
model and then treat the Cg distortion as a perturbation intermixing the trigonal CF terms. It 
was proposed in (see also [|ll|) and then confirmed in indirect experiment Q mentioned in 
the Introduction, that the ground state level and the first excited level form a pair of Kramers 
doublets 

\G±) = a|±l/2) ±6|^5/2) (27) 
\E±) = I ±3/2) (28) 

separated by small energy interval which was estimated as ~ 4.4 meV in a point charge approx- 
imation for the crystal field. Indeed, the purely electrostatic approximation is crude enough, 
because the covalent and exchange corrections should be taken into account as it is seen from 
our equations ( |2l| ) and ( [2^ ) . These corrections result in additional reduction of the CF splitting 
(see below), so one can be sure that the third CF level is high enough in energy (~ 14.0 meV 
according to calculations of |^), and the latter state is irrelevant to the low-energy excitation 
spectrum. 

To calculate the spectrum of RVB excitations one should diagonalize the mean field Hamil- 



tonian Hmf + SHmf given by equations (|l9| ) and (20) by means of Fourier transformation in 



Ce-sublattice which, in turn, has four sublattices in this orthorhombic structure with the use of 



the basis functions ( P7iP8| ). This is a cumbersome procedure in general case, but in the special 
case of CeNiSn one should take into account extreme anisotropy of magnetic response which is 
seen both in static magnetic susceptibility 0] and in quasi 2D character of neutron scattering 
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spectra |^. These experimental facts enabled us to presume that the magnetic anisotropy 
reflects the quasi 2D character of low-energy spin excitations, and, hence, to consider the phys- 
ical situation when these excitations possess the dispersion only in 6c-plane. The 2D character 
of spin-fermion excitations can be explained by the properties of indirect RKKY interaction 
which is responsible for spinon coupling. It is obvious that only those spins can be involved in 
the resonating valence bond which are connected by the antiferromagnetic exchange coupling. 
If one adopts the positive sign of the integral /pj^j^p in be plane {p is the distance between nn 
Ce ions), one can easily imagine that 15% lesser value of intersite distance in a direction [l^] is 
enough to have nearly zero or even negative value of oscillating function ^(kp, H-j — R-j+p) in the 
integral (0). One should emphasize that this assumption does not imply two-dimensionality 
of conduction electron spectrum which can hardly be imagined in CeNiSn lattice. 

Another possible source of anisotropy of magnetic excitation spectrum is the anisotropy of 
indirect exchange described by the factor BG{0,(f>) in eq. (|l^). The degree of anisotropy given 
by this mechanism can be estimated for a special case of the ground state ( p7| ) (see Appendix). 
This estimation shows that the "in plane" interaction exceeds the interattion along the a-axis 
when \a\ ^ \b\ in the case of small interionic distances and when \a\ <C \b\ in the case of 
great interionic spaces. Since the interionic distance in CeNiSn is neither big nor small and the 
experimental value of a can be estimated as o ~ 0.67 (see below), we have no real grounds for 
referring to this mechanism. 

One more reason for changing the sign of ^ is the anisotropy of the electron Fermi surface. 
However, the available information about the latter is rather scanty, so we have no firm basis 
for discussing this mechanism. 

In any case the assumption of two-dimensionality is not crucial for the results obtained, 
as soon as the exchange interaction along a-direction is not too big. The brief analysis of the 
influence of a-component of spinon dispersion on the neutron scattering spectra is presented in 
Subsection 5.4. 

Thus, assuming the 2D dispersion of spin liquid excitations, we project the Ce-sublattice 
onto be plane and find that the orthorhombic 2D elementary cell contains two Ce ions in the 
sites i = 1^ where ^ = 1,2 defines the sublattice (see figure 1). This network is defined by 
the Bravais vectors B = (6,0) and C = (0,c) and the basis vector d = (0, — 6/2,c/2 — O). 
Here O is the orthorhombic distortion which transforms one-ion hexagonal lattice into two-ion 
orthorhombic one. 

To describe the 2D spinon spectrum in a nearest neighbour approximation we define the 
coupling constants Tai in equation (|^) by two parameters, 

'^1 = = T\2,\'2 (29) 

for coupling within the same sublattice, and 

^2 = T\l^V2 = T\2,V1 (30) 



for intersublattice coupling. Similarly, two mixing integrals Bn' in equation ( |24| ) are introduced, 
i.e. the one-cite coupling constant Q 



Sgfi = B?A . (31) 



^This constant includes also the contribution of monoclinic distortion of electrostatic crystal field 
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and inter-site coupling constant 



Q2 



^GE 



11,1'2 



12,1'1 



(32) 



Such choice implies that the very existence of two sublattices is due to orthorhombic distortion 
O of trigonal lattice, so displacement of second sublattice (see figure 1) is, probably, the main 
source of the non-Coqblin hybridization. Then, introducing the Fourier transformation of spinon 
operator, 

= N-'/' E E E E ^l^'^i^^ k)/,5,r. exp (.kl) (33) 
1 ^=1 r 

and taking into account the fact that this type of displacement generates the crystal field which 
components are given by the Stevens operators (see also Subsection 5.2), we find the system 
of equations for the coefficients Hjj'''(^,k) and eigenvalues e^(k) , 

N 2 G,E ± 



EEEE^Prv.'(i^"' 

1' ^'=1 r' ly' 



(e', k) exp {ikl'} = e^{k)El'''{^, k) exp {ikl} 



(34) 



Here the coefficients -Dpp ^^/(R;;/) are defined as 



^IW(R = 0) = AcF, V. 



11 

GG;uu 
-.11 



(R = Rl)=Tl, V, 



GG;uiy(^ 
12 



Ri 



To. 



V'e'gM'^ = 0) = Gi V\,%,A^ = Ri) = 02 



(35) 



[see equations (29) - (^2[)], Ri is the coordinate of the nearest neighboring site. Then 
the secular matrix can be decoupled into two blocks, {{1G+),{2G+),{1E—),(2E—)} and 
{(IG—), (2G— ),(!£'+), (2£^+)} which correspond to two new Kramers doublets given by the 
following linear combinations, 



/k;.,±, = H^(l, k)/kl,G± + S;:(2, k)/k2,G± + S^(l, k)/kl,i^^ + H^(2, k)/k2,^ . 

Each block of the secular matrix has the form 

f TiMi(k) T2M2(k) Gi g2M2(k) \ 

TsMKk) riMi(k) gsMKk) Gi 

Gi G2M2{y) AcF y 

V G2M^{k) Gi Ace j 

The structure factors which define the spinon dispersion are determined as 

Mi(k) = 2 cos (kB) 

M2(k) = 4 exp 



(36) 



(37) 



kB-kC 



cos 



kB 



cos 



(¥) 



(38) 
(39) 



(B = (6,0) and C = (0,c)). 

The excitation spectrum is shown in figure 2 (the insert presents the dispersion of spin exci- 
tations in Coqblin-Cornut approximation (^1 = 0,G2 = 0) where the CF excitations are simply 
superimposed on the ground state spinon band). It is seen from this figure that the intermixing 
of spinon and CF excitations results in radical reconstruction of the low-energy part of this 
spectrum. We associate the characteristic temperature T* estimated from various experiments 
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as T* ~ 80K (see Introduction) with the characteristic energy scale of this structured part, so 
the values of 71 = 12. 7K and T2 = 18K were chosen to reproduce approximately the width 
of this energy interval. The values of other parameters {Qi, Q2 and Aqp) where taken to fit 
the positions of the main peaks in the neutron scattering spectra (see the next section). In 
particular, we adopted the value of 13.5 K for the renormalized CF splitting parameter A^^). 
Its reduction in comparison with the value of 48.4 K given by indirect experiment of [Q] is easily 
explained by the effects of covalent and exchange renormalization given by the terms B^, Bf 
and in eqs. (^Tj) and (p2[). It is obvious that the covalent repulsion B^ of the localized levels 
Ey from the free conduction band states above the Fermi level (2^,22) is governed by the value 



of hybridization matrix elements | V^y' I ^ • Since the Fermi surface of CeNiSn is dominated by the 
p-partial waves [^], the hybridization is largest for the jz = ±1/2 and jz = ±3/2 components 
of the f-states. Therefore, the ratio of hybridization matrix elements for |G±) ( p7[ ) and |-E±) 
(|28|) can be evaluated as | V^'q_|_P/| « < 1. This means that the negative covalent term 

reduces the CF splitting. The exchange contribution B^^ in (21) gives additional contribution 
to this reduction effect. 

The dispersionless behavior of spinon spectrum at the 6-facet of Brillouin zone [(0, 1/2, 0) < 
k < (0, 1/2, 1/2) ] is explained by the specific form of the spinon spectrum 

£^=1,2 (k) = 2 |ti cos (kB) ± 2T2 cos cos | (40) 

since this spectrum becomes dispersionless for kB = vr. 

It is important that inclusion of intermixing terms Qi^2 in the secular matrix ( |37[ ) does not 
alter the one-dimensionality of this branch (cf. upper and lower panels of figure 2) because 
the form factor M2 given by eq. (^) also turns into zero at kB = vr. Basing on the spinon 
dispersion law shown in this picture we will try in the next section to explain the neutron 
scattering data. The chemical potential of spin-fermions at T = which is determined by the 
constraint (^) is shown in this picture by the dashed line. 

Thus we find that the spin excitations in CeNiSn preserve to some extent the features of 
the one-site crystal field excitation spectrum: the multiband dispersion picture contains vast 
dispersionless parts which origin is the initial atomic crystal field splitting of the f-level. The 
density of spin-fermion states described by these bands is shown in upper and lower panel of 
figure 3. We see that instead of simple two- level picture of CF excitations there exists several 
peaks some of which can be considered as the remnants of CF levels but other reflect the van 
Hove singularities of 2D spinon spectrum. The spectrum demonstrates the pseudogap features 
presupposed in early phenomenological descriptions |^. However, it is remarkable that the 
chemical potential of spin-fermion excitations falls not into the hybridization gap but into the 
pseudogap between two van Hove peaks. 



4 Inelastic Neutron Scattering in Spin Liquid 

The inelastic neutron scattering provides unique experimental information concerning the struc- 
ture and dispersion of low-energy spin excitations, so these data look challenging for the general 
theory of spin liquid and especially for the case of interplay between HF anf CF excitations. To 
calculate the spectra of inelastic neutron scattering in the spin liquid, we start with the well 
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known equation for magnetic scattering cross section 



(fa 
dVLdE' 



2 \ ^ / 
76 \ K 



a/3 



(41) 



Here rrif, and e are the mass and the charge of electron, 7 = —1.91 is the gyromagnetic ratio 
for neutron, c is the light velocity. The wave vectors k and k' are momenta of incident and 
scattered neutron, respectively, and are the cartesian components of unit vector directed 
along the momentum transfer vector Q. The component 5q,/j(Q, Tloj) of the scatering function 
can be expressed as the sum over all possible initial (A) and final (A') states with the energies 
E\ and respectively. 



S«^(Q, fiu) = ^ (a I Q+ 1 A') (y \Qp\\)5 {hco + E^-E, 

AA' 



(42) 



{n\ is the probability distribution for initial state). Then, the operator Q of neutron-electron 
interaction can be expressed as a sum over all f ions. 



cell 

Q = ^ exp (iQRi) exp (iQdg) Qi^ 
1 € 

In the dipolar approximation the operator Q;^ turns into 



(43) 



(44) 



where g is the Lande splitting factor, -F(Q) is the ionic fomfactor and J is the total angular 



momentum operator. Thus the problem is reduced to calculation of the correlator ( A 



A 



In our model the source of neutron scattering are the transitions between the spin liquid 



excitations described by the operators (33) with the dispersion law e/x(k) given by the secular 
equation (|3^). Now the situation formally resembles the well-known paramagnetic scattering 
1 18 1 where the spin-fermion excitations A = k/x play the same part as the itinerant electrons 
in transition metals. Substituting the eigen vectors given by equation (|3^ ) in the scattering 
function (^2|), we find 



E^^' Ek kr'(k, Q)l * /r'(k, Q)n^(k) [1 - n^,(k - Q)] 6 {hco + e,{k) - e^{k - Q)) (45) 



where 



T'u' ( 



(e,k-Q) rz. 



(46) 



n^(k) is the Fermi distribution function for spinon excitations 

This expression can be simplified in case of CeNiSn. It is known from experiment [l^ that 
the dominant component which contributes significantly to the scattering function is 5aa(Q, hi-^) 
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component. Thus, for the momentum transfer Q in the b — c plain, which is the experimental 
condition for the most of experimental scans, one has for the scattering function 



^ ^ |/r'(k, Q)|' n^(k) [1 - n,ik - Q)] 6 {fiu + e^(k) - e^(k - Q)) (47) 

k 

Since the axis a is the axis of spin quantatisation, and since the operator Ja has only diagonal 
nonzero matrix elements in the basis of CF states, \Tv) (27,28), the expression for 11^^ (k, Q) 
has the form 



ir' (k, Q) = E E (^Qd^) {rF/ii, k)) * -^"^/{i, k - Q) (r. 

Tu i 



r 



Tv 



(48) 



Finally, taking into account the explicit equation ( pq ) for the spinon operator and expressions 
for the CF states (^^) one comes to the final equation for the matrix element /^^ (k, Q) 
(signs lb are for up and down partner of the Kramers doublet (p6|) respectively) 



/^^(k,Q) 



(s^(i, k)) * H^, (1, k - Q)eG± + (h^(i, k)) * H^, (1, k - Q)es^' 

(h^(2, k)) * H^, (2, k - Ct)@G± + (h^(2, k)) * H^, (2, k - CDQe^ 



+ 



X exp (iQd) 



where 



eci = ± 3a' 



®- = ^2 



(49) 

(50) 
(51) 



5 Inelastic neutron scattering spectra in CeNiSn 

In this section we compare the inelastic neutron scattering spectra calculated by means of eqs. 



(47) and (^) with dispersion law for spin-fermions given by secular matrix (^) and figure 2 
with the experimental data of the refs. 0, §, P, 10| 



5.1 Qualitative considerations 

The inelastic neutron scattering spectrum of CeNiSn in its gross features is formed by two main 
signals with the energies of ~ 2meV and ~ 4meV. It is seen from the lower panels of figures 2 and 
3 that these "peaks" can be ascribed to transitions from the occupied states in the lower band to 
the region of the empty levels where the states |G±) ( [27| ) and \E^) ( p^ ) are strongly hybridized. 
The most distinctly these features are seen in experiment for the momentum transfer vector 
Q in the 6c-plain and, hence, the magnetic response is connected with the x'"* component of 
dynamic magnetic susceptibility Q . As is known from the measurements of static succeptibility, 
the a-axis is the easy magnetization axis in CeNiSn |^]. Since the lower spinon band is formed 
mainly by the ground state doublet |Gib) ( p7[ ) where the components | it 5/2) play the main 
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part (see next subsection), the longitudinal susceptibility should be the strongest component of 
magnetic response at low energies, so the dominance of Q _L a momentum transfer vectors is 
naturally explained in our model. Then the most intense transitions ar expected for those final 
states where the greater fraction of | it 5/2) is admixed to the | it 3/2) states. 

The 2meV peaks are seen experimentally for Q = (0,0, 1) and Q = (0, 1,0) |, g. The 
corresponding momentum transfer vectors are indicated by solid arrows in figure 4 where the 
shaded area is the projection of Brillouin zone to the be plane). These processes correspond to 
the vertical transitions from the spinon "Fermi surface" to the region of strong hybridization 
near the point k = (0, 0, 1/2). The bold arrow marks these transitions in figure 2 (the integer 
reciprocal vector (0,1,0) or (0,0,1) are substracted). The same transitions are shown by the 
bold arrow in the lower panel of figure 3. The highest peak in DOS arises due to the van Hove 
singularities near the former CF excitation level. 

The peaks at E = 4meV are seen for Q = (0, 1/2 + n, Qc) (grey arrows in figure 4) where 
n are integer numbers |^ g. This maximum in the scattering cross section is due to the 
intraband transitions shown by grey arrows in figure 2. It is easily seen that since the 4meV 
peak corresponds to the transitions between two peaks of DOS (grey arrow in the lower panel 
of figure 3), and the 2meV processes are due to the transitions from the structurless part of 
DOS to the empty peak, the intensity of the former signal should be higher than that of the 
latter one. This result is also consistent with the experimental data ^. 

It worth to be noted that the general shape of the DOS in the lower band reminds the double- 
peak DOS with pseudogap postulated in early phenomenological models |^ |2|, ^. However, at 
least one more peak with higher energy should exist in the scattering spectrum given by the 
present model. This peak is formed by the interband transitions shown by the dashed lines 
in figures 2 and 3. Such peak with E ~ 7meV was found in the recent neutron scattering 
experiments [|lO| . 

Thus, the present model explains all qualitative features of the neutron scattering spectrum. 
To make the detailed quantitative description of experimental scans one should first specify the 
model parameters. 



5.2 Evaluation of the model parameters 

To calculate the spinon spectrum determined by the secular matrix ( ^7|) one should know the 
values of the intersite coupling constant Ti^2, the mixing parameters Qi^2 the energy Acf of 
CF splitting. To find the intensities of neutron scattering one needs also the magnitude of the 
orthorhombic distortion O and the values of parameters a, b which determine the structure of 
the ground state CF level |Gib) (|27|). All these quantities, except hybridization parameters 
Ghc and Qo are, indeed, predetermined by the experimental data available. The criterions for 
choosing the parameters 7i^2 and Acf are described in Section 3. The value of orthorhombic 



distortion O in units of lattice parameter c is 0.10 [17|. 

We have no independent experimental information to choose the values of mixing parameters 
Gi^2, so we use these parameters for fitting the experimental neutron scattering cross section. 
The best values, which gave us the possibility to describe the low-energy magnetic scattering 
function for all measured momentum transfer are Qi = —3.5 K, Q2 = —7.7 K. 

To estimate the parameters a, b entering the ground state CF level wave function (^) we 
used the experimental data on the anizotropy of low temperature static magnetic susceptibility 
at low temperatures P|. According to our model, this susceptibility should be mainly determined 
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by the Pauh-hke contribution of RVB excitations. 

The Pauh-hke contribution for magnetic field applied along the j-axis at low temperatures is 
proportional to the DOS of spinous N{sp) and to the square of effective (/j^^-factor of the spinous 



states |G±) at the Fermi level, Xaf^ ~ dx'' J\f{eF) {J\f{eF) is the density of spinon states at the 

Fermi level). The effective g^^^ -factor is determined as gj^^ = 2gj\{Gv\Jj\Gv')\ {v = it). Here 
gj = 6/7 is the Lande factor for Ce^"^ and Jj is the j- component of total momentum operator. 
According to (^) , the spinon states at the Fermi surface which determine the anisotropy of the 
susceptibility are described as the mixture of the bare wave functions |Gib) and \E±) (p^),(p8|). 
due to orthorhombic distortion of the trigonal lattice. As was mentioned above, we consider 
the orthorhombic distortion (see Fig. 1). Therefore, the wavefunctions in the vicinity of the 
Fermi surface can be approximated by 

\G±) = a^J[)\G±) + (jg^E^) (52) 

( P\ ( F) 

where Oav and Pav are the mixing parameters averaged over the Fermi surface. 

When evaluating the Pauli contribution, we neglect the anisotropy of spinon Fermi surface 
and assume that the ratio between the components of static susceptibility is given by the 
anisotropy of g factor: 

X" : : ~ bap : ISfeP : lOcf ■ 

( P) ( F) 

To evaluate the averaged mixing parameters aiiv and Pav we take into account that there 
are = 4 nearest neighbors from the second sublattice which take part in k-dependent hy- 
bridization. The averaging of the k-dependence of hybridization gives additional factor 1/2. 
The distance from the Fermi surface ep to the crystal field level is ~ 30 K according to our 
calculations (see figures 2,3). Therefore, the value of Pav at the Fermi surface can be estimated 

( F) 

as Pav ~ NQ /2{AcF — £f) j so Uav ~ 0.88. The similar value of parameter parameter aav can 
be extracted from the experimental data for the anisotropy of magnetic susceptibility 



2 



Xc V 3aa„a / 



(where si ~ 0.25 [^). One finds aav ~ 0.94 from ( |53| ) in agreement with the theoretical 
evaluation. The estimation of the magnetic susceptibility gives 



2 



^ = Xa ^ Xa ^ f agy (2 - 6a^) + 3 . 
° ~ Xc Xb \ 3a2^a2 j 



where 2 < sq < 3. Using (|54[ ) the value of |a| is determined as 



I 2ag, + 3 
3aL(2 + ^/so 



(55) 

so we obtain 0.65 < \a\ < 0.75. Since the contribution of |5/2)-component in the magnetic 
response is dominant even for \a\ ~ we actually deal with the | it 5/2) states in the neutron 
scattering experiments. It turns out that the neutron scattering cross sections weakly depend 
on the the value of a in the actual range of |o|, and we can not extract this parameter from the 
available experimental data. In the following calculations the value of \a\ = 0.67 is taken. 
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5.3 Quantitative description of constant-Q and constant-E neutron scatter- 
ing scans 

First, we see that our model reproduces quite well the "pseudogap" behavior of constant-Q 
scans first found in Figure 5 presents the theoretical curves for Q perpendicular to a axis 
in comparison with experimental data of ref. [@] and later results of ^ . It is easily seen from 
Figs 2 and 3 (lower panel) that this character of Q _L a scans is due to high density of states 
for the I lb 5/2) components in the region of strong hybridization around (0,0, 1/2) point. This 
psedogap-like behavior disappears for Q||a Q simply because the 2-meV peak is suppressed for 
this direction of momentum transfer in accordance with the geometrical factor in eq. ( pl| ) under 
condition Xaa > Xbb, Xcc- 

The next striking feature of experimental scans is the quasi ID character of 4-meV excitation: 
it was noticed ^ that the behavior of this peak for Q = (0,n + 1/2, Qc) practically does 
not depend on the value of Qc- Our next scan demonstrates, first, that this peak is indeed 
can be obtained in the spin-fermion model (figure 6, upper panel) and, second, that nothing 
is changed significantly in the shape of this scan with varying Qc (lower panel). The form of 
spinon dispersion curve explains this one-dimensionality: due to the dispersionless character 



of the energy band along the wave vector k in k = (0, 1/2, fee) direction [see equation (40) for 
explanation] , the contribution from the intraband transitions with E ~ 4meV to the scattering 
function (spanted arrows in fiture 2b) does not depend on Qc and thus mimics quasi one- 
dimensional scattering. 

One more remarkable property of experimental constant-Q scans is their modification along 
the (0, Qb, 0) direction. It was found (see |9|, figure 3) that the form of these scans change from 
the structure with definite 4.2 meV peak at Q^ = 1.5 to that with 2 Mev peak at Qi, = 1. 
It is seen from figure 7 that our theoretical scans reproduce this variation quite satisfactorily, 
including the featureless shape of the scan at intermediate value of Qh = 1.2. 

The theory explains also why the 2 meV peak which is distinctly seen for the vertical 
transitions (0, Qf,,0) and (0,0, Qf,) (see figures 4,5) practically disappears for the "diagonal" 
direction Q = (0, 1, 1) (figure 16 of |p). This kind of anisotropy is explained by the fact 
that the phase factor exp (iQd) in (^) is close to unity for Q = (0, 1, 1) due to the small 
orthorhombic shift O = 0.10 (see Subsection 5.2). Then, since the condition Qg± ~ ©Eq: is 
valid for the chose value of a = 0.67, the matrix element (^8[), (49|) practically turns out into the 
orthogonality relation between the different quantum states (3^) with the same reduced wave 
vector. This is why the intensity of scan at Q = (0, 1, 1) is by order of magnitude smaller than 
that of Q = (0, 0, 1) or Q = (0, 1, 0). 

Finally, we discuss the characteristic features of constant-E scans. The calculated scans for 
E = 4.2meV are presented in figure 8 together with the experimental data taken from |^]. It is 
no wonder that the theoretical scans are in satisfactory quantitative agreement with experiment: 
it simply follows from above explanation that the maxima of (0,(5fe,0) and (0,(5fe, 1) scans at 
Qb = 0.5, 1.5 are connected with the fact that the 4meV peak arises just at Q = (0, n -|- 1/2, x) 
with integer n and any x. 



More puzzling feature of the constant-E scans was noticed in [10| where it was observed 
that the peak at Qh = 1.5 for the (0, (5b,l) scans with E = 4.2meV splits into two peaks 
for lower energy transfer E = 3.3meV and E = 2.5meV. This feature is reproduced, at least 
quantitatively, in our calculations (see figure 9). Both the appearance of the doubled peak and 
the increase of the splitting for lower energies is seen in the theoretical curves. Even the higher 
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intensity of the left peak is reproduced in the calculated scan. To investigate the reason of this 
asymmetry we calculated also the constant-Q scans with Q = (0, Qh, 1) for Qt, = 1.25, 1.5 1.75. 
Our figure 10 demonstrates that for higher energy transfer the intensity of Qb = 1.5 response is 
higher than that of Qf, = 1.25, 1.75 whereas for lower energy transfer the situation is reversed. 

Thus, the totality of calculated scans shows that not only the general shape but also the 
anisotropy of spinon spectrum correlates with the details of experimental neutron scattering 
scans. These results demonstrate clearly that the anisotropic magnetic response of CeNiSn is 
due to anisotropy of the spectrum of low-energy spin-fermion excitations in this system. 

5.4 Influence of 3-dimensionality of the spinon spectrum 

To check the robustness of our results to assumption about mainly two-dimensional character of 
spinon dispersion we recalculated the spinon DOS and inelastic neutron scattering intensities. 

We introduced in the secular matrix for the spinon spectrum two more parameters, i.e. the 
"overlap integral" = Tu y^ = 7i2,i'4 where the indices (3,4) stand for the Ce ions above 
and below the sites (1,2) respectively, in a-direction, and Q3 = B^l^,^ = B^^,^ with the same 
notations. These parameters result in further doubling of the crystal cell lattice and further 
removing the degeneracy of spinon spectrum and appearance of spinon dispersion in a-direction 
of the Brillouin zone. As to the spinon DOS, the overlap parameter is responsible mainly 
for for broadening of the second van Hove peak (-25 K peak in the lower panel of fig. 3) and 
the influence of the hybridization parameter ^3 is seen in broadening and lowering of the 3 
K peak in the mentioned DOS. But the general shape of the DOS, and, therefore, the form 
of the inelastic neutron scattering spectra for Q in 6c-plane are practically not modified at 
least for 73 and G3 less then « AK, i.e. more then 50 per cent of the characteristic values of in- 
plane interaction parameters. Correspondingly, the neutron scattering peaks for the momentum 
transfer in be plane are also slightly broadened, but the general quantitative agreement is still 
satisfactory, although the CF splitting parameter should be slightly changed. If one chooses, 
e.g., T3 = = 3.5K, the renormalized value of Acf = 19. 5K should be taken to get an 
agreement with the experimental data for all scans demonstrated in figures 5-10. 

However, the introduction of third component in spinon spectrum provided us one more 
possibility to check the model assumptions, i.e. to calculate the neutron scattering for Q out of 
6c-plane. In fig. 11 the constant-E scan for Q = (Q^, 3/2,0) is presented in comparison with the 
experimental data of [^. The oscillating character of scattering intensity with Qa is distinctly 
seen in experiment and is shurely reproduced in the calculated curve. It is interesting that the 
shape of the theoretical curve is more complicated than the rough sinusoidal guideline shown 
in fig. lid of Ref. and the attentive look on the experimental points confirms this picture at 
least qualitatively. 

When calculating the theoretical curve we weighted the inelastic cross section only with 
the form factor F{Q) [see eq. (^5[)]. This procedure resulted in better agreement with the 
experiment than that including also the orientational factor 1 — Q'^. This means that the 
contribution of Imx*^ and Im^"^"^ to this cross section is also significant due to, e.g., the zigzag 
distortion of Ce chains an a-direction. 
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6 Low-temperature specific heat of CeNiSn 



To be sure that our approach gives really universal description of low-energy spin excitations in 
CeNiSn, we calculated the magnetic contribuion to specific heat of CeNiSn within the approxi- 
mate mean-field approach to RVB thermodynamics with the use of the same parameters for the 
spin fermion spectrum which were chosen for description of the neutron scattering data. Since 
the samlpes studied in the neutron scattering measurements, apparently, were not as perfect as 
the recently grown samples which demonstrate the metallic type conductivity and low value of 
7 = C/T ~ 40mJ/mole (5[, we used for comparison the data for the sample #2 from Q. The 
value of 7(T = 0.1 K) = 80mJ/mole was chosen to eliminate the ultralow-T upturn due to 
nuclear and impurity contribution to the magnetic specific heat. The results of our calculations 
are presented in fig. 12 (solid curve). The experimental data presented by diamonds in figure 
12 are obtained from those for that sample by cancelling the phonon contribution to specific 
heat. The latter was considered to be the same as in LaNiSn crystal ||20t| . 

When comparing the experimental and theoretical curves for the specific heat, one should 
take into account the contribution of conduction electrons and add it to the magnetic specific 
heat given by the spin fermion excitations. We have taken it to be linear in T within the low-T 
interval < T < 10 K with the Sommerfeld coefficient 7e=18 mJ/mol K^. This means that 
the electron contribution to the heavy fermion density of states in good samples can be as large 
as ~ 50 per cent in agreement with the theoretical predictions for the spin liquid state of the 
Kondo lattices |22]. It is seen from figure 12, that the theoretical curve describes quite well the 
temperature behavior of C/T = j + bT at low temperture and gives the maximum of C/T at 
T = 6 K. The latter is observed for the more perfect samples |2^ where the data for higher 
temperatures are available. 

It should be noted that the contribution of spinon dispersion in o-direction (see Section 5.4) 
practically does not change the form of C(T) curve since the van Hove peaks survive in 3D 
spectrum, and the spinon chemical potential is still in the dip of the DOS. 

Our previous semi-phenomenological studies have shown that the model of heavy fermion 
liquid with the pseudogap in the spin-fermion excitations induced by the interplay with the CF 
excitations describes quantitatively the low-temperature specific heat and thermal expansion of 
CeNiSn IQ, 11]. Now we see that the consequent microscopic theory gives the detailed picture 
of the low-energy excitations and low-temperature thermodynamics of this system which is in 
reasonable quantitative agreement with the experiment. 



7 Concluding remarks 

To conclude, we have found that the general theory of spin-fermions acquires some specific fea- 
tures in the low-symmetry Kondo lattice materials with soft CF excitations, such as CeNiSn and 
CeRhSb. These CF excitations modify the spectrum of spin excitations in a very characteristic 
manner, and as a result, the density of spin-fermion states has several peaks due to interaction 
between the spin-fermion excitations and CF excitations. Because of existence of these peaks 
the spin response of CeNiSn and CeRhSb reminds that of spin-gap materials. 

Whereas the largest peak in spinon DOS at ~ 2 meV above the spinon Fermi level appears as 
a result of interplay between the HF and CF excitations, all other singularities in one-spinon and 
two-spinon DOS arise due to quasi two-dimensionality of spin excitation spectrum (although 
these singularities are robust enough against the appearance of the third component of spinon 
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dispersion). The source of this anisotropy of the excitation spectrum is the anisotropy of RKKY 
interaction in a real space. It is worth mentioning that the standard slave-boson description 
used in a previous phenomenological model ||2|, |l^ also predicts highly anisotropic spectrum of 
HF excitations due to the sf-hybridization anisotropy in k-space. 

If the levels | =F 5/2) dominate in the ground state Kramers doublet |Gib) (27) (see Sec- 
tion 5.2), then the effective intersite overlap integral T-^'^ is determined in this theory by the 

combination of hybridization matrix elements, ^A^A'/(^i^ ~ ^f)i which turns into zero for 
k = [1.0.0] (see, e.g., p3|). This means that the f-electron wave functions form the non-bonding 
states along this direction. Moreover, the expression ^^'0*^0/^^ ~ ^f) which determines in 
this case the mixing integrals B^^^ turns into zero for k = [1,0,0] due to the similar reasons. 
Thus one can expect that the main characteristic features of the interplay between the HF 
and CF states, including the pseudogap, would develope in the 6c-plane of the Brilloun zone. 
However, this kind of theory has an undesirable consequence: in the mean-field slave-boson 
approximation the spin and charge pseudogaps simply coinside, and this statement obviously 
contradicts the experimental data for CeNiSn. 

Another version of the mean-field description of pseudogap in CeNiSn was presented recently 
in [^]. It was postulated in this theory that the ground state of Ce ion is the doublet | =F 3/2), 
and the CF excitations where not taken into account. The pseudogap in the HF spectrum 
appeared due to the mentioned non-bonding character of f-f overlap matrix element T'^^(k). 
Then the low-temperature behavior of specific heat and magnetic susceptibility are determined 
by this pseudogap, and the satisfactory agreement between the theoretical description of these 
quantities and the experimental data for CeiSn can be achieved. 

As to the inelastic neutron scattering cross section, the theoretical results presented in [24] 
are too scanty to conclude about agreement between the theory and the experiment. One can 
definitely mention the qualitative agreement, at least for the constant-E scans (figure 11 of 
Ref . |24(] which is more or less similar to our figure 8). As to the constant-Q scans, the slave- 
boson description is able to describe the pseudogap in the two-particle HF DOS (figure 10 of 
Ref. [p^ ) , although the experimental picture is, apparently, much more reach then the smooth 
curves presented in that figure. One should mention also, that the theoretical spectrum for 
Q = [1/2,0,0] presented in figure 10b demonstrates the peak in a-direction which is absent 
in experiment because of the geometrical factor entering the neutron scattering cross section 
for this direction [|^, It is unclear also whether this theory is able to describe the quasi 
ID character of neutron scattering along the [0,n -|- 1/2, Qc] direction. On the face of it, this 
task seems to be rather difficult since the matrix element which governs the HF dispersion law 
behaves as V'^{k1 + ky){\ + 15A;^) (see equation 2.10b in Ref.p4[). Besides, the CF level scheme 
adopted in this model is not supported by the experimental data available [^. Finally, the 
coincidence of the charge and spin pseudogap in the mean field slave boson theory mentioned 
above does not allow to reproduce the metallic charater of low-temperature resistivity. The 



absence of T in the theoretical curve p{T) is clearly seen from figure 13 of Ref. [24 



The inelastic neutron scattering method is an adequate tool for observing the spin-liquid 
excitations in the heavy-fermion materials. These measurements give the detailed information 
concerning the density of states of spinon particle-hole excitation spectrum, and the dispersion 
of various branches of spinon excitations can be restored by analysing various constant-Q and 
constant-E' scans. The general picture of low-energy excitation spectrum obtained from this 
analysis allowed us to describe also the low-temperature specific heat of CeNiSn (see also 
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All other important characteristics of this material, such as magnetic susceptibility, thermal 
expansion, NMR relaxation rate, etc, also can be described quantitatively within the same 
approach. 

The specific features of spin-liquid excitation spectrum do not violate the conduction electron 
spectrum which remains essentially metallic, although the anisotropy of spin excitations can 
result in anisotropy of scattering mechanism for conduction electrons, and the soft crystal field 

excitations can influence their paramagnetic response. In any case, CeNiSn and related materials 
should be treated as metallic Kondo lattices with specific type of spin excitations. 
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A Appendix 



To calculate the anisotropy factor in Cornut-Coqblin Hamiltonian one should expand the Bloch 
waves entering the RKKY interaction (pT|),(p!3D in partial waves entering the hybridization 
integrals in equation (^) and then superimpose the direction connecting two ions with the 
quantization axis by means of the matrices of finite rotations. 

In the case of short enough interionic distances one can neglect the influence of angular 
dependence of the partial waves on the values of radial overlap integrals. This approximation 
]T^ . in a simplest case of the doublet formed by the j = 5/2 states with the same jz value, 
i.e., {Ti/) = \j zizM) results in the following equation for the factor of anisotropy of exchange 
interaction, 

Bg{9) = {BM{e)f , 



BM{e)= (-iri2/ + i)^rp/(cos0) 

«=0,2,4,6 

Zt' = diO{l, M + 1/2) + rilO{l, M - 1/2) 



= V(7 + 2M)/14, ,7M = V(7-2M)/14, 0(/, j) = Cfo^Cg^ (56) 

Here C^Mim Clebsch-Gordan coefficients. This equation can be found in [0] for the case 

oi 6 = 0. Since the CF wave functions transform in accordance with the spherical group the 
only source of anizotropy is the definite direction of quantatization axis which lowers symmetry 
to the axial one. Therefore the interaction depends on the polar angle 6 only. 

More complicate is the anisotropy factor for a mixed representations of the type of eq. (27). 



In this case the wave functions transform according to the irreducible representation of the 
trigonal group with third order rotational axis parallel to the quantization axis. In this case 
the RKKY infraction becomes ^-dependent, 

Bcie, (/)) = 1 1 + - (5 - 9a2) P2(cos e) + -{l + a^) Pi{cos 9) 

"-P|(cos0)cos(3(/>) I (57) 



7V10 

For the specific cases of interaction "in plane" {Q = vr/2) and along the quantization axis 
{6 = 0) the factor B does not depend on the angle (p. Therefore the ratio between the "in 
plane" interaction and interaction along z-axis can be expressed in terms of single parameter a 
entering the wave function (p7|), 



BG{e = 7r/2) J 1 25 + 81a2 



BG{e = Q) [40 4-3a2 



(58) 



It is seen from (^) that Bg{9 = tt /2) / Bg{0 = 0) ^ 1 only when |l/2)-compomemt of the 
crystal field wave function ( p7|) is dominant. 

In the case of the large interionic distances the asymptotic behavior of RKKY exchange 



integrals is dominated by = ±1/2 partial waves |26], and the "in plane" interaction is stronger 
than that along the quantization axis {9 = 0) in the case when the |5/2)-component of the crystal 
field wave function (1271) is dominant. 
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Figure captions 



Figure 1: The be plain of CeNiSn lattice. Two Ce sublattices are denoted by black and grey 
circles, resp. The orthorhombic distortion O (solid arrows) transforms the simple hexag- 
onal lattice into two-sublattice orthorhombic one. The vector d is the basis vector of 
two-ion elementary cell 

Figure 2: The spinon dispersion in CeNiSn along the high symmetry directions of 2D Brillouin 



zone obtained for the following paramenter values (in Kelvins) of the secular matrix (37): 
7i = 18.0,72 = 12.7, AcF = 13.5, = 3.5,^2 = 7.7. Horizontal line at the energy 
£' = — 13 K is the Fermi level position at T = 0. The arrows indicate the intraband and 
interband transitions responsible for the main peaks in neutron scattering cross section 
(see text). 

Insert: the non- hybridized spinon spectrum with CF level (dased line) superimposed 

Figure 3: Upper panel: the density of states of unhybridized spinon spectrum 

Lower panel: the density of states for spin-fermions hybridized with CF excitation. All 
other designations are the same as in fig. 2 

Figure 4: The Q vectors corresponding to the main transitions in spinon bands in 2D reciprocal 
lattice. The shaded areas are the projection of BZ onto be plain and the regions in BZ 
where the final states are located 

Figure 5: Comparison of experimantal (diamonds) and theoretical (curves) constant-Q scans. 
Upper panel: Q=(0,0,1.2). Experimental poins are taken from 0. 
Lower panel: Q=(0,1,0). Experimental points are taken from ^] 

Figure 6: Theoretical (lines) and experimental (points) constant-Q scans with 4meV peak. 
Upper panel: Constant-Q scan for Q = (0,1.5,1). Experimental points are taken from 
lH) Hi (squares) and [|lO| (diamonds). 

Lower panel: Constant-Q scans for Q = (0,1.5,1) (solid line), Q = (0,1.5,0.75) (dotted 
line) and Q = (0, 1.5, 0.5) (dashed line) 

Figure 7: Constant-Q scans for: Q=(0,1.5,0); (0,1.4,0); (0,1.3,0); (0,1.2,0); (0,1.1,0); (0,1,0) 
(curves 1 -6, respectively). Position of 4meV and 2meV peaks are indicated by arrows. 
The vertical scales are shifted for clarity 

Figure 8: Theoretical (curves) and experimental (dots) constant-E scans (E = 4.2meV) along 
Q = (0,(5f),0) (squares and solid line) and Q = (0,Qfe,l) (diamonds and dashed line) 
directions. Experimental points are taken from 

Figure 9: Theoretical (curves) and experimental (dots) constant-E scans along (0, Qb, 1) direc- 
tion: (a) E=4.2mV; (b) E=3.3meV; (c) E=2.5meV. Experimental points are taken from 

lie 



Figure 10: Constant-Q scans for: Q = (0,1.25,1) (dashed line); Q = (0,1.5,1) (solid line); 
Q = (0,1.75,1) (dotted line). 

Figure 11: Theoretical (curves) and experimental (dots) constant-E scans {E = 4.25meV) 
along Q = (Qa, 3/2,0) direction. Experimental points are taken from M 
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;ure 12: Theoretical curve for low-T contribution to the specific heat in comparison with 
the experimental data taken from Q (diamonds). The constant electronic contribution 
7e; = 18 mJ/mol is added to the theoretical curve (see text) 
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